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Abstract Variation in age structure and body size 
benefits are identified to understand the evolution of life 
history. Here, we estimated the age structure and body 
size of two species of odorous frogs (Odorrana margaretae 
and Odorrana grahami) by using skeletochronology. The 
ages at sexual maturity of O. grahami and O. margaretae 
in both sexes were 1 and 2 years, respectively. For both 
sexes, the maximum age observed in O. margaretae was 
six years. For O. grahami, the maximum age observed in 
males and females were 4 and 5 years, respectively. Males 
and females did not differ in mean age in the two species. 
The average body size of both species considerably 
differed between sexes, with females being larger than 
males. The body size of females was also larger than that 
of males when the effect of age was removed. We also 
found positive correlations between body size and age 
within each sex in O. margaretae, but only for female in 
O. grahami. The female-biased sexual size dimorphism 
of the two species suggested that fecundity selection for 
larger female size may increase the reproductive output. 
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1. Introduction 


The global population declines of amphibians have highlighted 
the lack of knowledge in population dynamics (Green, 2003). 
The population dynamics are underpinned by mortality 
rates which are also affected by life-history traits (eg, growth 
rates, age at sexual maturity and longevity; Wells, 2007). The 
knowledge of individual lifespan and population age structure 
facilitates the identified factors that affect the mortality 
of certain life-history stages (Biek et al., 2002) and develop 
appropriate conservation strategies for anurans (Wells, 2007). 
Hence, identifying the patterns of age structure may help us to 
understand life-history traits in organisms (Zhong et al., 2018). 

The sexual size dimorphism (SSD), which is the difference 
in body size between males and females, is widespread (Wells, 
2007; Adams et al, 2020), and females are the larger sex in 90% 
of anuran species (Shine, 1979). Previous studies indicated that 
the body size of individuals differs significantly between sexes, 
and the body size of females is larger than that of males in most 
anurans (Leclair and Laurin, 1996; Liao and Lu, 2010; Liao et al., 
2011; Liao and Chen, 2012). This pattern has been successfully 
output favours large body size in females, thereby leading to 
female-biased SSD (De Lisle and Rowe, 2013; Liao, 2013; Liao et 
al., 2015). 

Skeletochronology represents a reliable method for 
determining age structure in anurans on the basis of presence 
of the lines of arrested growth (LAGs) recorded in the cross- 
sections of phalanges (Castanet and Smirina, 1990). Over the 
past 30 years, this method has been used in estimating the age 
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of anurans from different areas (temperate: Hemelaar and 
Van Gelder, 1980; Guarino and Erismis, 2008; Kyriakopoulou- 
Sklavounou et al., 2008; subtropical: Morrison et al., 2004; Liao 
and Lu, 2011; Liao et al, 2011; Zhang et al., 2021; tropical: Khonsue 
et al, 2000; Kumbar and Pancharatna, 2001) without killing 
frogs. 

In the present study, we investigate the age structure 
and body size of two species of odorous frogs (i.e, Odorrana 
margaretae and Odorrana grahami) inhabiting the ravine streams 
in China by using skeletochronology. The green odorous frog 
(O. margaretae) is a species distributed in Hengduan Mountains 
where it is restricted to a narrow range of elevations, ie, 500- 
1650 m (Fei and Ye, 2001). The disk-fingered odorous frog (O. 
grahami) is one of the most widely distributed species in China 
with elevations elevation ranging from 1720 to 3200 m (Fei and 
Ye, 2001). The two species have similarities in habitat types, diet 
and activity time (Fei et al., 2010). However, their population 
size has declined in recent years due to human activity. 
Although altitudinal variations in body size of O. grahami have 
been reported by Li et al. (2013), information about the age 
structure data of the two species during the breeding period is 
unavailable. Our aims are to compare age structure, body size, 
and the relationship between age and body size between sexes 
in the two species, where these demographic parameters are 
poorly known, and provide conservation suggestion on the 


basis of the age structure and body size of two species. 
2. Materials and Methods 


2.1. Study area The population of O. margaretae was located 
in a stream of Wawushan National Park in Ya’an City 
(29.64°N, 103.09°E), at an altitude of 1105 m a.s. in the Hengduan 
Mountains of Western China. The population of O. grahami 
was located in a stream of Midu County in Dali City (25.12°N, 
100.61°E) at an altitude of 1550 m a.s.l. in the Hengduan 
Mountains of Western China. The two study areas have a 
strong seasonal variation in climate. The monthly average 
temperature ranges from 3 °C to 22 °C in Wawushan National 
Park and from 11 °C to 22.5 °C in Midu County. 


2.2. Sampling methods All frogs were captured in ravine 
streams from late June to early August 2018 during their 
reproduction. We collected 43 O. margaretae individuals in 
June 2018 and 43 O. grahami individuals in August 2018. We 
determined their sexes and sexual maturity by observing the 
secondary sexual characteristics (ie, nuptial pads in males and 
evidently swollen belly due to egg masses in females). We used 
a calliper with accuracy up to the nearest 0.1 mm to measure 
body size (snout-vent length, [SVL] within each individual. We 
then clipped the third phalanx of the longest toe of the right 
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hind limb, and stored them in 4% neutral buffered formalin for 
skeletochronology. 


2.3. Skeletochronology Following the methods used for 
determining the age of each frog (Liao and Lu, 2010), we 
removed the skin and muscle tissues of each digit and used 
water to wash the remaining bones for 2 h. The digits were 
decalcified in 5% nitric acid for 48 h, washed in running tap 
water overnight and stained with Ehrlich’s haematoxylin for 
100 min. Subsequently, we dehydrated the stained phalanges 
of each digit through successive ethanol stages of 70%, 80%, 
95%, and 100% for approximately 1 h at each concentration. 
Phalanges were processed for paraffin embedding in small 
blocks (50 mm x 15 mm x 10 mm; L x W x H). We used a 
rotary microtome to obtain thin cross-sections (13 ym) of the 
mid-shaft diaphysis, selected the phalanx with the smallest 
medullar cavity, and then mounted them onto glass slides. We 
used a light microscope to observe the sections and the Motic 
BA300 digital camera mounted on the Moticam2006 light 
microscope at 400 x magnifications to photograph the best 
section. LAGs were identified by two persons (M. J. Zhong and 
L. J. Yuan), who were equally experienced in the technique and 
identification criteria and agreed on the final age estimation. 

The number of LAGs corresponded to the age of individuals 
due to distinct temperature cycles experienced through the 
year (Morrison et al., 2004). We considered the surface of the 
bone as an additional LAG based the specimens collected 
during the breeding seasons after emergence from hibernation 
(Liao and Lu, 2010). We used the presence of the Kastschenko 
Line (KL, the interface between the endosteal and periosteal 
zones; Rozenblut and Ogielska, 2005) to confirm the endosteal 
resorption of LAGs for the study species by following the 
protocol by Castanet and Smirina (1990). 


2.4. Statistical analysis We performed all statistical tests by 
using the SPSS 22.0 software (Statistical Product and Service 
Solutions Company, Chicago). We used parametric tests 
when the data met assumptions of the tests (ie. normality of 
the distribution and homogeneity of the variance) and used 
non-parametric tests when assumptions were not met. We 
first identified the age structure between sexes by using the 
Chi-square test and then used the Mann-Whitney U-test 
to compare differences in body size and age between males 
and females. We conducted ANCOVAs and treated SVL as 
a dependent variable, and age as covariate to test significant 
difference in body size between sexes when correcting the 
age effects. Finally, we used general linear models and treated 
SVL as a dependent variable, and age and sex as fixed factors 
to assess the effect of interaction between age and sex on SVL. 
All values given were shown as mean + SD, and the level of 
significance was P < 0,05. 
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3. Results 


3.1. Skeletochronology A total of 80 individuals were aged 
successfully by skeletochronology. We obtained 32 males and 
9 females in O. margaretae and 18 males and 21 females in O. 
grahami (Table 1). Six individuals lacked LAGs were excluded 
from the analysis. Extremely closely spaced hematoxylinophilic 
lines (so-called double LAGs) were found in two O. margaretae 
females. We observed false LAGs (ie, LAGs that were faint and 
did not form a complete ring in the bone section) in 1 male and 
1 female in O. grahami, However, we did not find problems in 
estimating their age. Moreover, we found that the endosteal 
resorption of the innermost LAG occurred in 7.5% of males and 
8.5% of females in two species based on the occurrence of KLs 
(Figure 1). 


3.2. Age at sexual maturity and longevity The age at 
sexual maturity was estimated to be two years in males and 
females in O. margaretae. For O. grahami, the minimal age of 
males and females at sexual maturity was 1 year. The lifespan 
of individuals in O. margaretae was estimated to be at least six 


years in males and females, whereas the lifespan of males and 
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females in O. grahami did not exceed 4 and 5 years, respectively 
(Figure 2). 

Males and females did not differ in age in O. margaretae 
(Mann-Whitney U-test: Z = 0.469, n, = 32, n, = 9, P = 0.639) and 
O. grahami (Z = 0.320, n; = 18, n, = 21, P = 0.749). In O. margaretae, 
the average age of adult males and females were 3.88 + 1.21 


(range = 2-6) and 3.90 + 1.24 (range = 2-6) years, respectively. In 
O. grahami, the average age of males and females were 2.57 + 140 
(range = 1-4) and 2.39 + 0.98 (range = 1-5) years, respectively. 


3.3. Age structure The adult age of O. margaretae ranged 
from 2 years to 6 years in males and females (Figure 2). Age 
distributions differed significantly between sexes (Chi-square 
test: y” = 12.902, P < 0.001). The adult age of O. grahami ranged 
from 1 year to 4 years in males and from 1 year to 5 years in 
females (Figure 2). Age distributions did not differ significantly 
between sexes (y” = 0.231, P = 0.631). The predominant age class 
of males and females was 4 years (34.19%) in O. margaretae and 3 
years (30.8%) in O. grahami. 


3.4. Body size The average SVL differed significantly between 
males and females in O. margaretae (Mann-Whitney U-test: Z = 
23.116, n, = 32, n, = 9, P < 0.001) and O. grahami (Z = 2.338, n, = 18, 


Table 1 Body size and age of two species of odorous frogs (Odorrana margaretae and Odorrana grahami) in southwestern China. Values in 


descending order are mean + SD and sample size. 


Odorrana margaretae 


Odorrana grahami 


Sex 


SVL (mm) Age (years) SVL (mm) Age (years) 

Female 87.87+16.90 3.90+1.24 72.24+19.68 2390.98 
n=9 n=9 n=21 n=21 

Male 75.54+6.19 3.884+1.21 5951225 2.57+1.40 
n= 32 n= 32 n=18 n=18 


Figure 1 Two hematoxylin-stained cross sections of the phalangeal bone (A: a 3-yrs old male Odorrana margaretae; B: a 4-yrs old 
female Odorrana grahami) in southwestern China. Arrows indicate the arrested growth lines (LAG). KL represents Kastschenko Line, the 
interface between the endosteal and periosteal zones. Scale bar: 100 pm. 
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n, = 21, P = 0.019). When correcting the age effect, the differences 
in body size between sexes remained significant in O. margaretae 
(ANCOVA: F, 4, = 220.676, P < 0,001) and O. grahami (F, x = 
6.275, P < 0.019). A significant effect of interaction between age 
and sex on body size in O. margaretae (F, 4 = 763.382, P < 0.001) 
revealed that the relationship between age and size differed 
between males and females. For O. grahami, no significant effect 
of interaction between age and sex on body size (F, 4 = 1369, P 
= 0.250) was observed, showing that the relationship between 
age and size did not differ between sexes. 

The average SVL (Mann-Whitney U-test: females, Z = 
10.868, n, = 9, n, = 21, P < 0.001; males, Z = 18.901, n, = 32, n, = 18, P 
< 0.001) and average age (females, Z = 16.430, n, = 9, n, = 21, P < 
0.001; males, Z = 21.519, n, = 32, n, = 18, P < 0.001) of O. margaretae 
females and males were larger than those of O. grahami females 
and males. 


3.5. Relationship between age and body size Regression 
equations revealed significant relationships between SVL and 
age within each sex in O. margaretae (females: y = 11.303x + 
43.823, F, o= 7.534, r, = 0.518, P = 0.029; males: y = 2.693x + 64.554, 
F, 3, = 9.425, r, = 0.239, P = 0.005) and in female O. grahami (y 


12 
10 


Number of individuals 


Age (yrs) 


oF N WRN QA Yo 


Figure 2 Adult age structure (male: close bars; female: open 
bars) of (A) Odorrana margaretae and (B) Odorrana grahami in 
southwestern China. 
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= 9.003x + 49.084, r, = 0.249, F, », = 13.181, P = 0.002), but not in 
male O. grahami (y = 4.180x + 49.528, r, = 0.097, F, 1s = 2.008, P = 
0.176). A large overlap of SVL amongst the estimated age classes 
occurred within each sex, showing that age cannot be predicted 
from body size (Figure 3). 


4. Discussion 


We investigate for age structure and body size in O. margaretae 
and O. grahami. We provide clear evidence that females have 
larger body size than males when correcting for the age effect, 
but no significant difference is observed in the average age 
between males and females in the two species. We also find 
that body size displays a positive correlation with age within 
each sex in two species except male O. grahami. Our results 
suggest that the relatively short lifespan in both species is likely 
to result in increasing mortality rates, thereby decreasing the 


population size. 
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Figure 3 The relationship between age and body size of (A) 
Odorrana margaretae and (B) Odorrana grahami in southwestern 
China (female: close circle, full line; male: open circle, dotted 
line). 
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The use of phalangeal growth marks under low 
environmental temperature during the hibernating period 
depresses physiological functions of the frogs and can assess 
individual age in most anurans (Morrison et al, 2003; Lai et al., 
2005; Huang et al., 2014; Yu et al., 2018). The formation of LAGs 
in O. margaretae is affected by inactivity for a long period (Li et 
al, 2013). In this study, the hibernation in the two species begins 
in mid-October, and the inactive period lasts until early May, 
which produces the clear LAG formation for all frogs. The 
accuracy of the age determination in the frogs is often affected 
by double and false lines (Hemelaar and van Gelder, 1980). Here, 
we find 2 individuals with double lines and 2 individuals with 
false lines which do not affect age estimation and data analysis. 

As critical transitions in life-history traits, the age at sexual 
maturity and longevity remarkably display the influences of 
individuals’ fitness (Charlesworth, 1980; Xiong et al., 2020). For 
O. margaretae, the age at sexual maturity of males and females 
arise after the second hibernation due to the low temperature 
experienced at high altitude. Similar to some anurans (Esteban 
et al, 1996; Guarino et al., 2003; Liao and Lu, 2010; Liao and Lu, 
2012), the age at sexual maturity in O. grahami in both sexes 
was | year, suggesting that individuals possibly spend time 
in reproduction due to high temperature at low altitude. We 
find that the maximum ages recorded are 6 years in males and 
females for O. margaretae, 4 years in O. grahami males, and 5 
years in O. grahami females. Relatively short longevity for both 
species suggests that high mortality rates are likely to linked to 
human activity in studied populations. 

SSD may be consequences of differences in most key life- 
history traits between sexes (e.g, age structure and growth rate 
before the attainment of sexual maturity; Monnet and Cherry, 
2002; Zhang and Lu, 2013). Sex differences in age at maturity 
do not seem to explain differences in SSD in the studied species 
due to the same age at maturity in males and females. For Rana 
limnocharis, the non-significant interaction between sex and 
age on body size shows that age-SVL relationships between 
sexes do not differ in slope, and males and females have similar 
growth patterns, which do not affect SSD (Liao et al., 2011). We 
find that SSD in O. grahami is not affected by growth rate due 
to the non-significant effect of interaction between age and sex 
on body size. By contrast, the growth rate in O. margaretae is a 
major factor affecting SSD because the interaction between age 
and sex on body size differs significantly between sexes. 

The SSD is regarded as an outcome of sexual and natural 
selection (Fairbairn, 1997). We find that independent from age, 
females are larger than males in the two species, suggesting 
that natural selection on larger body size in females favours 
the improved fecundity or offspring size (Trivers, 1974). Sexual 
selection can drive the evolution of SSD, in which increased 


male body size may promote male mating success (Andersson, 
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1994). However, we do not provide the evidence that male body 
size is positively associated with male mating success in the 
two species. Finally, differences in mortality rates (e.g. due to 
predation) between sexes may lead to differences in the average 
age of adults and may affect variation in SSD (Liao et al., 2015). 

Body size changes with altitude and/or latitude in most 
anurans, with individuals living in high-altitude and/or latitude 
having larger body size than those living in low altitude and/ 
or latitude (Miaud et al., 1999; Lu et al., 2006; Matthews and 
Miaud, 2007; Chen and Lu, 2011). A previous study on body 
size variation across four O. grahami populations has suggested 
that in disagreement with Bergmann’s rule, body size does not 
change along the 830 m altitudinal gradients (Li et al., 2013). We 
do not find increases in body size in the 1550-m population with 
altitude when we analyze data on life history across five O. 
grahami populations. 

In most anurans, age displays positive correlation with 
body size within each sex (Esteban et al., 1996; Ryser, 1996; 
Kyriakopoulou-Sklavounou et al., 2008; Liao and Lu, 2011). 
We find that body size is positively correlated with age for 
both sexes in O. margaretae. However, in other anurans a 
positive correlation may be true for only one sex (Gibbons 
and McCarthy, 1984; Leclair and Castanet, 1987; Cherry and 
Francillon, 1992). For O. grahami, SVL displays a significant 
correlation with age in females but not in males. 

In combination with life-history traits, data on age structure 
in anurans may provide important insights into demography 
within and amongst populations, which have an important 
implication for the decisions of conservting species. The age 
structure collected on O. margaretae and O. grahami is being 
used to develop stage-structured population models to estimate 
population dynamics and test for the causes of population 
decline in the two species. 
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